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ABSTRACT 

Deviations from the gravitational inverse-square law would imprint scale-dependent features 
on the power spectrum of mass density fluctuations. We model such deviations as a Yukawa- 
like contribution to the gravitational potential and discuss the growth function in a mixed 
dark matter model with adiabatic initial conditions. Evolution of perturbations is considered 
in general non-fiat cosmological models with a cosmological constant, and an analytical ap- 
proximation for the growth function is provided. The coupling between baryons and cold dark 
matter across recombination is negligibly affected by modified gravity physics if the proper 
cutoff length of the long-range Yukawa-like force is > 10 h~ 1 Mpc. Enhancement of gravity 
affects the subsequent evolution, boosting large-scale power in a way that resembles the effect 
of a lower matter density. This phenomenon is almost perfectly degenerate in power-spectrum 
shape with the effect of a background of massive neutrinos. Back-reaction on density growth 
from a modified cosmic expansion rate should however also affect the normalization of the 
power spectrum, with a shape distortion similar to the case of a non-modified background. 
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1 INTRODUCTION 

General relativity has passed many importan t tests up to the lengt h 
and time scales of the observable universe iPeeblesll2002L Eo04) . 
The body of evidence is quite impressive, considering the enor- 
mous extrapolation from the empirical basis, but the gravitational 
inverse-square law and its relativistic generalization are supported 
by high-precision tests from measurements in the laboratory, the 
solar system and the binary pulsar only up to scales of < 10 13 cm 
jAldel bereer et all2003l) . 

Despite the overall success of general relativity, non-standard 
theories have received much recent attention, largely motivated by 
finding alternative explanations for the 'dark' sector of the uni- 
verse. On one hand, MOND proposals make gravity stronger on 
scales of galaxies to explain flat rotation curves without dark mat- 
ter jMilgrom|["l983l : ISanderslI 19981) : on the other hand, alternative 
theories aim to find a mechanism for cosmic acceleration without 
dark energy, as a result of gravity leaking into higher dimensions on 
scales comparable to the horizon iDvali et al]|2000l) . Such radical 
modifications of gravitational physics may seem difficult to test, but 
our ability to make precise observations of cosmological dynamics 
on scales of > 10 /i _1 Mpc is now good enough to permit detailed 
tests of gravity theories even on these scales. 

Examples of such work include Iwhite & Kochanekl 1200 ll) . 
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who determined constraints on the long-range properties of grav- 
ity from weak gravitational lensing cosmic shear. Constraints on 
Newton's constant from t he prim ordial abundances of light ele- 
ments were discussed by lUmezu et alJ 120051) . The implications 
of MOND on large scale structure in a Friedmann-Lemaitre- 
Rober tson- Walker (FLRW) universe were discussed in iNusseJ 
(2002). Recently, proposals have been made to extend these the- 
ories to a general covariant framework l Bek ensteinl l2004) and 
evolut ion of pert urbations in this theory has been considered 
ISkordis et all2005l) . Finally, large galaxy surveys have been anal- 
ysed looking for signatures in th e power spectrum due to d evia- 
tions from the inverse-square law JSealfon et alj2005tlShirata et alJ 
2005). In the present paper, we shall concentrate on this last aspect. 



We follow a number of recent authors, who have considered 
additional contributions to the gravitational poten t ial in the form 
of Yukawa-like terms fwhite & Kochanekl 1200 ll : IPeeblesI l2002t 
ISealfon et alll2005t IShirata et alj|2005h . We consider the growth 
function and the power spectrum of matter density fluctuations in 
these Yukawa-like models; we extend previous work by removing 
a restriction to fixed cosmological parameters and we also explore 
some possible back-reactions from a modified cosmic expansion 
rate. We still adopt the cold dark matter (CDM) model for structure 
formation with adiabatic initial conditions, but we also consider 
variants due to hot dark matter. The paper is structured as follows. 
In Section|2| we introduce Yukawa-like contributions to the gravi- 
tational potential and study the linear perturbation equation and the 
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Figure 2. Linear growth rate, normalized to the pure Newtonian case Da 
(Act = 0), as a function of D\. The background cosmology is a flat 
ACDM model with r2M0 = 0.25, whereas the Yukawa-model parame- 
ters are fixed to Act = 0.5 (i.e pAa ~ 0.27) and A = 10 Mpc. The lines 
from above to below correspond to k = 10~ e , 10 — 3 , 10 — 2 , 10 — 1 and 
1 Mpc , respectively. The growth function evolves as Da at very early 
times or small scales, and as D.^ PAa when a > Xk. 



Figure 1. Linear growth rate as a function of the scale factor a in a flat 
ACDM model with Qmo = 0.25 for Yukawa-model parameters Act = 0.5 
and A = 10 Mpc. Extra growth is effective at late times and on large scales. 
Upper panel: dependence on scale. Extra growth is activated for a > kX. 
The full, short-dashed and dashed lines are for k = 10" 2 ,5xl0- 2 and 
10 Mpc , respectively. Lower panel: fractional residuals between the 
above numerical results and the fitting formula. 



growth function. In Section[3] we give some insights on the coupled 
growth of baryons and CDM. In Section|4] we illustrate the effect 
of modified gravity on the power spectrum, whereas Section [5] is 
devoted to the interplay between mixed dark matter and gravity 
physics. Growth of perturbations in a model with modified cosmic 
expansion is addressed in Section|6| SectionQcontains some final 
considerations. 



2 LINEAR PERTURBATION EQUATION 

The dynamics of linear perturbations can be easily extracted from a 
Newtonian approach by writing down the fundamental equations of 
fluid motion (e.g. Peacock 1999, chapter 15). Usually, acceleration 
is related to mass density via Poisson's equation, but here we allow 
for deviations. We model the weak field limit of the gravitational 
potential, cf>, as a sum of a Newtonian pote ntial and a Yukawa-like 
contribution (see Whit e & Kochanekl200lL and references therein): 



(1) 



The factors ctNwt and ayi™ account for deviations of the effec- 
tive gravitational coupling constant from Gn, the value actually 
observed on small scales. </>Nwt is the solution of Poisson's equa- 
tion, whereas ^kw, which describes forces mediated by massive 
particles, is the solution of the modified Poisson equation; in co- 
moving coordinates, this can be written as 

V 2 (^Ykw </>Ykw . „ r, ,s 

2 rj- = 47rG N /9M<5(x,t), (2) 

with a the cosmic scale factor, pm the total matter density, 5 the 
mass density contrast and A a proper length cutoff. The potential in 
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Figure 3. Fractional deviation of the growth function, at the present time 
a = 1, between an Einstein-de Sitter (EdS) model and a flat ACDM model 
with SImo = 0.25. The Yukawa-model parameters are Act = 0.5 and 
A = 20 Mpc. The suppression due to a cosmological constant is scale 
dependent. 



Eq. Q is oc 1/x both on small scales (a; <C A/a), with an effective 
coupling constant G = (aNwt + ctYkw)GN, and on large scales, 
with G oo — ctNwt Gn - Putting 



CtNwt 



1 + Aq 



(3) 



and ctYkw = -Aq, <j) is in numerical agreement with results on 
the scale of the solar system, whereas Goo = (1 + Aq)Gn. For 
Aq > 0(< 0), gravity is enhanced (suppressed) on large scales. 
The potential in Eq. can be derived i n a relativistic gravity model 
that obeys the equivalence principle iZhvtnikov & Nesteilll994l) . 
Although the Yukawa-like form is quite specific, all long-range de- 
viations can be characterised empirically by an amplitude, Act, and 
a length scale, A. 

The basic equation for the linear evolution of the Fourier com- 
ponent 5k of the density fluctuation, with comoving wave-number 
k, is 



5 k + 2H(t)& = 



(4) 
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where dot denotes time derivative and H = a /a is the time de- 
pendent Hu bble term. In Fourier space, the p otential in Eq. Q is 
jwhite & KochanekfcOOltlSealfon et all2005h 



Ok = 75 y\G>l d\ ,' QNwt 3 QYkw ) , 



with 

^ = QfNwt + 

= 1 + Aa 



1 + {a/a x y 
(a/a\) 2 



(5) 



(6) 



(7) 



1 + (a/a A ) 

and a\ = Xk. Fr om Eqs. J4| [5), we obtain feealfon et aT]|2005t 
IShirata et alj2005l) 



o k + 2HS k + 



3 ^ImoHq 



y{a/ay, a?Nwt, «Ykw)5k = 0, 



(8) 



where Qmo is the present density of total matter in units of the 
critical one and Ho is the Hubble constant, parameterized in the 
following as 100 h km s _1 Mpc _1 . A deviation from the inverse- 
square law thus causes the rate of growth of the Fourier amplitudes 
to depend on wavelength. The inclusion of a Yukawa-like contribu- 
tion imprints the scale A/ao on the power spectrum. Let us examine 
the growth between an initial epoch, aj n , and the actual time, an, in 
more detail. On a very large scale, i.e. k <C am/A and y ~ 1+ Aa, 
modes always experience the modified growth due to the effective 
gravitational constant. Intermediate scales, initially below the max- 
imal comoving cutoff length A/a m , will cross the scale for modi- 



fied growth at some time ai n 



ax 



ao. Finally, very small scales 



k 3> ao A will evolve in a Newtonian potential. As a result, the 
growth of fluctuations is no longer independent of scale, even at 
low redshift. Some examples of the growth function are plotted in 
Fig-d 

Equation 0, with <j) given as in Eq. has been solved an- 
alytically for an Einstein-de Sitter model (E dS, SImo = L^ko = 
0) in terms of hypergeometric functions in lShirata et alj (2005 J. 1 
Here, we want to examine approximate solutions in the general 
case with cosmological constant. Let us begin with the EdS model. 
A Yukawa-like contribution with Aa > 0(< 0) enhances (sup- 
presses) the growth of fluctuations on scales sufficiently larger than 
A. The growth law becomes 



PAa 



V25 + 24A7 



(9) 



The growth function, D, of the CDM+baryons density fluctuations 
can be approximated by smoothly interpolating between the ana- 
lytic results. We propose, for the growing mode, 

pao/ c i 



D 



Ykw(a; Aa, A) = &n— \ 1 + 



a/a in 



1 + C2ax/a in 



(10) 



where the normalization has been chosen to be D — 5; n a/ai n at 
early times a ~ a,i n and the coefficients ci and C2 represent a fit to 
numerical evolution. For ai n = 10~ 3 , we get 



ci = 1.94088 
c 2 = 1.42350 



0.170668Aa 



(11) 
(12) 



The 

10~ 2 



fitting formula 

< k/(h Mpc" 1 ) 



works quite well, for 
< 1, for the parameter range 



1 We refer to SectionltHfor considerations about modified cosmic expansion 
rate and adopt here a standard FLRW framework. 



,1+PAo 



-0.25 < Aa < 1 and 2 < A/(/i _1 Mpc) < 100, with fractional 
residuals always under 1%. 

The growth function changes at late times when Qyio 7^ 1 
since, after matter ceases to dominate the expansion rate, fluctua- 
tion growth halts. This must be considered together with the late- 
time effect due to deviations from the inverse- square law since the 
comoving scale for modified growth, A/a, moves to smaller values 
with increasing time. It is easy to account for both of these effects. 
Let Da denote the growth function for a model with non-zero cur- 
vature and cosmological constant but in the absence of Yukawa-like 
terms. Then, by the replacement of a/ai n with Da in Eq. 1 101 . we 
can approximate all the dependence of the scale-dependent growth 
function on time, curvature and cosmological constant in presence 
of a modified gravitational potentia l. A u seful approximation for 
Da can be found in ICarroll et alJ ll992l) . In the lower panel of 
Fig.Q we compare the numerical evaluation to the fitting formula. 
The accuracy is quite high, with performance better than 1% for 
a parameter space similar to the EdS case. The fitting formula ap- 
proximates the asymptotic limit on very large scales or late times, 
a\, as D A +P&a . This is a natural generalization of the 
form that applies for the EdS model. Although the result 
is no longer exact if Qao 7^ 0, it is very accurate in practice. In 
Fig. [2] we can see how the linear growth index changes from the 
usual Newtonian trend to an enhanced index for decreasing wave- 
numbers and decreasing redshifts and how well the £)^ +PA ° ap- 
proximation works at late times. 

In the presence of a Yukawa-like term, the suppression due to 
a cosmological constant with respect to the EdS model is most ef- 
fective at large scales, in contrast to the Newtonian case, where the 
suppression is scale-independent. Let us consider the case Aa > 0. 
Whereas a Yukawa-like contribution is effective at late times and on 
large scales, extra growth is inhibited by a cosmological constant 
when matter becomes sub-dominant in the dynamics. In Fig.|3] we 
plot the fractional deviation in the growth function between an EdS 
cosmology and a flat ACDM model. 

The hypotheses used above are standard in modelling the ef- 
fect of a Yukawa-like contribution to the gravitational potential on 
the growth function ISealfon et all2005l:ISnirata et alj2005T) . In the 
next sections, we will relax some assumptions and extend the anal- 
ysis to mixed dark matter models. 



3 BARYONS 

Baryons imprint characteristic features into the power spectrum 
as the direct result of small density fluctuations in the early uni- 
verse prior to recombination, when they are tightly coupled with 
the photons and share in the same pressure-induced oscillations 
lEisen stein & HvJl998l) . As far as the Yukawa-like potential affects 
all kinds of matter in the same way, the modified coupled equations 
for perturbation growth in a mixture of baryons and CDM can be 
written as 



5b 
5c 



+ 2H 



6c 



(13) 



4nG^Ly 

i IMO 



- k 2 /(kjy) 

^B0 



^co 



where Sb and 8c are the density fluctuations in baryons and CDM, 
respectively, fcj = a^/farGrTpM / c s is the comoving Jeans wave- 
number in a Newtonian potential and c s is the sound speed. On a 
small scale, baryon pressure causes the growth rates in baryons and 
CDM to differ. Gravity-driven growth is not allowed for baryons 
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above a critical wave-number, if k/Qtjs/y) 2> 1. This condition 
can be re-written for small values of Aa as 



fe » fa 1 + 



Aa/2 



1 + (fcjA/a) 2 



+ C(Ac 



(14) 



A positive value of Aa pushes the domain of predominance of 
pressure over gravity to smaller lengths. When pressure over- 
comes gravity, a Yukawa-like term does not affect the dynamics 
of baryons. As usual, their main behaviour will be oscillatory, with 
slowly declining sound waves described by the WKB solution (e.g. 
iPeacockl 19991 equation (15.60)). 

Either on large scales or after the Compton drag epoch, when 
the sound speed drops by a large factor, pressure effects are negligi- 
ble. Solutions of Eq. 1 1 3i are simple since the matrix in the driving 
term has time independent eigenvectors: (1, 1) and (fico, — Obo)- 
Hence, a coupled perturbation will be quickly dominated by the 
fastest-growing mode with Sb = 8c and baryons quickly catch up 
with the CDM. 

Since we have been considering a proper length 
A > 10 h~ Mpc, the interaction of matter and radiation at 
decoupling is not affected by modified gravity: the cutoff at this 
time is > 10 /i _1 Gpc in comoving units, so its effects appear 
on a scale too large to be observed. Prior to the drag epoch, the 
comoving baryonic Jeans scale, Aj = 2n/kj, tracks the sound 
horizon, which is of the order of 100 /i~ 1 Mpc at a ~ 10 -3 . 
Hence, for a proper cutoff length A > 10 /i _1 Mpc, the effects of 
a Yukawa-like term on the effective Jeans wave-number, Eq. 1 1 4-1 . 
are really negligible. After the drag epoch, the oscillatory mode 
of baryons matches onto a mixture of pressure-free growing and 
decaying modes. Such modes are given by the solutions of the 
perturbation equation discussed in Section|2|for a Qyio = 1 matter 
dominated universe. Because the form of the WKB solution is 
preserved and because the modes after recombination are affected 
by modified gravity only at later times, velocity overshoot is still 
produced with the output after the drag epoch being maximised 
when the wave consists of a pure velocity perturbation. 

Although the position of the acoustic peaks is nearly un- 
affected, later evolution can still affect the scales of baryonic 
oscillations for small values of A and large values of Aa. A 
Yukawa-like term can alter the height of peaks, with the first 
peak enhanced with respect to the second one when Aa > 0. 
Baryonic aco ustic oscillations have bee n proposed a s probes of 
dark energy <Blake & Glazebrookll2003l : ISeo & Eisensteinll2003t 
lMatsubara l2004). but most of the attention is on position rather 
than height, so that the effect of modified gravity is negligible at 
the pr esent status of accuracy of galaxy redshift surveys ICole et all 



4 POWER SPECTRUM 

The linear power spectrum can be constructed from the transfer 
function as 



P(k,a) oc k n T 2 (k) 



D(z, a) 

D(Z, din) 



(15) 



where T(k) is the matter transfer function at the drag epoch, and n 
is the initial power spectrum index, equal to 1 for a scale invariant 
spectrum. As discussed in the previous section, we can assume that 
standard theory provides the transfer function at the last scattering 
surface. In what follows, we wil l use the fitting formulae for T(k) 
in lEisenstein & Hij <1998tll999h . 
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Figure 4. Effect of a Yukawa-like term on the power spectrum. Upper panel: 
the full, dashed and long-dashed lines are for Aa = 0,0.3 and 0.6 re- 
spectively, with A being fixed at 10 h~ x Mpc in a flat ACDM model with 
h = 0.7, SImo = 0.25, Qbo = 0.05 and n = 1. Units of the power 
spectrum are arbitrary. An enhancement of gravity (Aa > 0) boosts large- 
scale power. Lower panel: fractional residuals between the pure Newtonian 
case of the upper panel and a re-normalized power spectrum in a flat model 
with h = 0.7, !1mo = 0.30, Hbo = 0.05, n = 1 and Yukawa-model 
parameters Aa = 0.19, A = 5 h~ 1 Mpc. The effect of enhanced gravity 
is similar to that of lower matter density. 



A positive long-range contribution to the gravitational poten- 
tial enhances clustering on scales larger than A. The shape distor- 
tion of the power spectrum can be quite significant: see Fig. |4] A 
Yukawa-like contribution with Aa > changes the spectrum in 
a way that resembles a lower matter density, with a near degener- 
acy between the pair {Aa, A} and Qmo- An illustrative example 
is contained in the lower panel of Fig. [4] Larger fractional residu- 
als appear at the scale of the acoustic peaks, since in our example 
we have varied SImo but kept Q,bo fixed. As a matter of fact, the 
location of baryonic oscillations is nearly independent of Aa and 
could break the compensation. As well as a change in density pa- 
rameters, a Yukawa-like contribution can mimic a variation in the 
primaeval spectral index, with values of n < 1 nearly degenerate 
with Aa > 0, over a large range of wave-numbers. 

The degeneracy between matter density and Yukawa parame- 
ters could be quite significant when determining cosmological pa- 
rameters, since power- spectrum analyses provide one of the most 
important direct constraints on the value of £Imo, whereas other 
methods that are nearly independent of a Yukawa-like contribution, 
i.e. methods based on the cosmic microwave background radiation 
or supernovae of type la, are mostly sensitive to other combinations 
of cosmological parameters. 



5 MIXED DARK MATTER 

The introduction of massive neutrinos or other forms of hot dark 
matter couples the spatial and temporal behaviour of perturbation 
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Figure 5. Scale-dependent growth function in a MDM cosmology, eval- 
uated at the present epoch. We consider a flat ACDM model with h = 
0.7, S!mo = 0.25. Full lines represent Newtonian cases (Act = 0); thick 
lines include neutrinos (3 species and f v = 0.10); short-dashed lines ac- 
count for a Yukawa-like term with Act = 0.20 and A = 12 fc _1 Mpc. 
Scale dependent features due to massive neutrinos are nearly degenerate 
with modified gravity. The long-dashed curve stands for the growth func- 
tion with Act = —0.20, A = 12 h _1 Mpc, 3 neutrino species and 
fv = 0.05. Effects due to neutrinos can be cancelled by a Yukawa-like 
term with Act < 
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Figure 6. The effect of massive neutrinos on the shape of the power spec- 
trum is nearly degenerate to that of a Yukawa-like contribution. We consider 
a flat ACDM model with h = 0.7, f2 M o = 0.25, Q B o/^MO = 0.15 and 
n = 1. Upper panel: the full line represents the power spectrum of the den- 
sity weighted CDM, baryon and neutrino perturbations for f v = 0.05, 3 
neutrino species and Newtonian potential (Act = 0); the dashed line stands 
for Act = 0.19 and A = 10 /i Mpc without neutrinos. The units for 
the power spectrum are arbitrary. Lower panel: fractional residuals after re- 
scaling. 



growth jM<Jll996llHu & Eisenstei3ll998l:lmsenstein & HuJll999l) . 
On scales smaller than the free-streaming length, neutrinos do not 
trace CDM and baryon perturbations, which in turn grow more 
slowly because of the reduction of the gravitational source. As neu- 
trinos slow down and their free-streaming scale shrinks to even 
smaller scales with increasing time, the transfer function acquires a 
non-trivial time dependence. 

As has been discussed in iHu & Eisensteinl ll998l) : 
lEisenstein & HvJ Il999t) . the transfer function in a mixed dark mat- 
ter (MDM) cosmology with CDM, baryons and massive neutrinos 
can be decomposed into a scale-dependent growth function that 
incorporates all post-recombination effects, and a time independent 
function that represents the perturbations around recombination. 
The late-time evolution of density-wei ghted CDM and baryo ns 
perturbations can be accurately fitted by iHu & Eisensteii] 1998 ) 



D cb (y,k) = <1 + 



Q , r \ Pcb/0.7 



D A (V) 



l-Pcb 



(16) 



where Da is the growth function in absence of neutrino free- 



streaming, y 



'i I a 



with a cq the epoch of matter-radiation 



equality, yt s is the free-streaming epoch as a function of scale, 
and p c b = (5 — \/l + 24/ c b) /4 with / c b the matter fraction in 
CDM+baryons. The fitting formula in Eq. I16i interpolates across 
yf E between the suppressed logarithmic growth rate at small scales, 
which is 1 — p c b in a EdS model, and the usual growth rate 
at sufficiently large scales. An analogous formula for the growth 
rate for the CDM+baryon+neutrino perturbations can be found in 
lHu&Eisensteir] | 1998). If we replace Da with the growth function 
-D Y kw> discussed in Sec|2| D c k will contain all the dependence of 
the transfer function on time, curvature, cosmological constant and 
Yukawa-like contribution in presence of neutrino free-streaming 2 . 
An amplitude Aa > has a similar effect to that of a non-zero neu- 
trino mass, both reinforcing large-scale evolution. Some examples 
of growth functions are plotted in Fig. [5] Scale-dependent features 
due to neutrinos can be almost perfectly cancelled by a Yukawa- 
like term. 

It turns out that features in the power spectrum due to the 
transition between free-streaming and infall of massive neutrinos 
have much in common with imprints induced by enhance d grav- 
ity: se e Fig. [6] where the fitting formulae from lEisenstein & HvJ 
( 1999), which do not account for baryonic oscillations, have been 
used. An enhancement of gravity (Aa > 0) can distort the shape 
of the power spectrum in nearly the same way as a non-zero neu- 
trino mass, both boosting large-scale power. Even at low redshifts, 
some scales can be still either in the free-streaming regime or be- 
low the A-cutoff, so there is an interplay between the two effects. 
A true deviation from the inverse-square law can be interpreted as 
a significant neutrino fraction /„ of the total matter density: see the 
lower panel of Fig. [6] In a low density flat model with three nearly 
degenerate neutrino species and /„ < 0.10, an adequate approxi- 
mation for such a degeneracy, is, if A < 50 ft -1 Mpc, 



(/») 



0.312 - 0.252 



(Aa) 



(17) 



with the proper length A in units of h^ 1 Mpc. For a smaller cutoff 
proper length, a larger deviation Aa is needed to mimic the effect 
of massive neutrinos. This compensation works up to the scales 



Note that the growth functions in IHu & EisenstemHl998t) have been nor- 
malized to D = a/a e q at early time, so that Dyi<w must be conveniently 
re-scaled. 
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Figure 7. The growth function, evaluated at the present epoch, in a model 
with modified background expansion. The shape of scale-dependent fea- 
tures due to a Yukawa-term are negligibly affected by a modified cos- 
mic rate. Upper panel: the full line stands for a modified EdS model with 
Aa = 0.5 and A = 20 Mpc; the short-dashed line represents the growth 
function in a non-modified EdS model with Aa = 0.4 and A = 21 Mpc. 
Lower panel: fractional residuals after re-scaling. 



where nonlinear or redshift-space effects begin to a ffect the shape 
of the power spectrum, i.e. for k < 0.1 h Mpc -1 l Meik sin et alJ 
Il999l) . 

The fractional contribution of neutrinos to the total mass den- 
sity has been derived from the 2dFGRS data lElsar0v et al. 2002). 
Adding prior information from independent cosmological probes, 
lElgaretv et all <2002l) found /„ < 0.13 (at 95% confi dence), i.e an 
upper bound on the total neutrino mass < 1.8 eV. Hspergel et all 
(2003) combined 2dFGRS data with measurements from WMAP 
and obtained a tight limit m v < 0.69 eV. These results are str ongly 
dependent on the assumed priors l Elgar0v & Lahav 2003J) and 
could be significantly affected by a long-range deviation from the 
inverse-square law. Since a long-range enhancement could equally 
lead to extra large-scale power, the upper bound on the neutrino 
mass would be over-(under)estimated if Aa > 0(< 0). 

This degeneracy applies only at one epoch, and the evolu- 
tion of P(k) will be different in the two cases. Although neutrinos 
mainly affect the growth function, they also contribute a residual 
suppression of power in the transfer function on small scales at the 
drag epoch. On the other hand, a Yukawa-like contribution, which 
is characterized by a proper length cutoff, is only effective at later 
times. 



6 MODIFIED COSMIC EXPANSION 

In the previous sections we have followed the usual approach 
by assuming that a modified gravitational potential only af- 
fects fluctuations while leaving the background cosmology in- 
tact <White & Kochaneljl200lk INusse3l2002t ISealfon et allEoO^: 
IShirata et alJ 120051) . We now want to address the effect of a 



perturbed cosmological expansion rate, due to modified grav- 
ity physics, on the dynamics of large-scale structure. In fact, 
Lagrangian-based theories of modified Newtonian gravity, in which 
the FLRW background cosmology still holds in the absence of fluc- 
tuations, have been proposed l Sanders 2001 ), but such frameworks 
are still controversial. Possible cosmological implications of new 
long-range gravitational interactions have been discussed within 
the framework of Newtonian cosmology. D'Olivo & Rvan 
showed that for a homogeneous, isotropic universe, with a Yukawa- 
type term in the gravitational interaction, a Newtonian cosmology 
exists and the equations of motions are the same as for the ordi- 
nary Newtonian case, i.e. the Friedmann equations. The only dif- 
ference is that the effective gravitational constant at large distance, 
Goo, enters the equations instead of the local measured value, Gn. 
Although no proof has been produced that such predictions from 
a Newtonian cosmology can be used in full relativistic cosmology, 
this modified FLRW model with an effective coupling constant pro- 
vides a motivated framework where implications of a modified cos- 
mic expansion rate on density fluctuations can be explored. 

We should therefore consider using the standard form of 
Friedmann equations, with Goo in place of Gn, and the standard 
cosmological parameters, with the only change being to normal- 
ize energy densities to a re-scaled cosmological critical density, 
= 3Hq/ (87rGoo). Hence, the usual cosmological param- 
eters are related to the new ones by 



Gn 

Goc 



(*Nwt 



(18) 



Within this model, local estimates of the matter density parame- 
ter, as those from cluster galaxy mass and luminosity, should be 
re-scaled when being compared to values determined from cosmic 
distances. 

A Yukawa-type potential causes a change in the cosmic ex- 
pansion rate of the background. The effect on linear perturbations 
is easily understood. Now, also the Hubble drag term in Eq. de- 
pends on the effective coupling constant Goo, whereas the form of 
the driving term on the right hand side is not affected. A positive 
long-range interaction (Aa > 0) enhances the drag term in Eq. J4j, 
so that growth is suppressed with respect to the non modified cos- 
mological dynamics. We rewrite the linear perturbation equation 
as 



5 k + 2HS* + 



with 



0. 



HNwt 



(19) 



H(a) = Ho V"mo W « 3 + "Io w + "ko w « 2 - (2°) 

As can be seen after comparison with Eq. {SJ, Eq. < 191 can not 
be simply obtained by substituting each cosmological density pa- 
rameter Qio with the corresponding f2^ kw , since the driving term 
is suppressed by a factor l/awwt due to the new role played by 
the effective coupling constant in the background dynamics. Let 
us consider the equivalent of an EdS model, i.e a flat model with 
= 1. The perturbation equation can now be solved in term 



of hypergeometric functions, 

-3/2+ Pfj 



5k = G] 



(-) 



(21) 
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+ 



2F1 



Vp 5 pp 1 / a \ 
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Tp/3 5 pp 9 / a \ 2 

X ' Fl [a'4 "f^"'''^) _ < 

where 2-F1 is the hypergeometric function, Ci and C2 are normal- 
ization constants determined by the initial conditions and pp = 
(5 - 725 _ 24/3) /4 with /3 ee Aa/(1 + Aa). The second 
term in Eq. <2H is the growing mode which interpolates between 
5 oc a 1 "* 1 '* at small scales (a <C oa) and the usual S oc a at 
late times and/or large scales. Due to the modified background, a 
Yukawa-like contribution with Aa > 0(< 0) determines suppres- 
sion (extra-growth) at small scales, whereas there is no variation on 
large scales. 

The global effect on the shape of growth function is actually 
similar to the case of a non-modified cosmic expansion rate dis- 
cussed in the previous sections, the main difference being an over- 
all normalization factor. The two cases are really indistinguishable 
when the differences in the growth index between large and small 
scale, i.e. pac and pp, are equal. Let us further explore the near- 
degeneracy between modified and non-modified cosmic expansion 
by considering the EdS model in the two cases. We recall that the 
different definitions of critical density, ^mo™ = 1 or ^mo = 1 cor- 
respond to different values of pyio- A value Aa in a non-modified 
background distorts the power spectrum in nearly the same way 
as Aa' = (3Aa - Ap\ a )/ [±p\ a + 3(1 - Aa)] in a modified 
cosmology with approximately the same cutoff proper length. In 
Fig.Q we show the degeneracy in the shape of the growth func- 
tion between modified and non-modified cosmic expansion rate. 
Differences show themselves only in the normalization factor, with 
growth of fluctuation suppressed for enhanced long-range gravita- 
tional interaction, due to the boosted cosmic expansion. 



7 CONCLUSIONS 

We have discussed the dynamics of large-scale structure in the pres- 
ence of a Yukawa-like contribution to the gravitational potential. 
The strength parameter and its cutoff length can introduce pecu- 
liar scale-dependent features in the shape of the power spectrum of 
density fluctuations. The main imprint of enhanced gravity is en- 
hanced large-scale power, in a way that resembles the effect of a 
lower matter density. If the extra long-range force couples to all 
kinds of matter in the same way, then baryons and CDM evolve 
together after recombination. Furthermore, for A > 10 /i _1 Mpc, 
the transfer function at the Compton drag epoch is not significantly 
affected by Yukawa-like contributions. 

With regard to evolution, enhanced gravity acts very similarly 
to hot dark matter: both effects yield a scale dependence of the 
growth function, even at low redshifts, with a suppression at small 
scales. Leakage of gravity on a proper scale of A > 10 /i~ 1 Mpc 
could nearly completely erase the signatures of non-zero neutrino 
mass on the shape of power spectrum, making the detection harder 
and softening the upper limit on the mass. 

A modified cosmic expansion rate should not affect the scale- 
dependent distortion of the power spectrum shape due to Yukawa- 
like contributions. Whereas the exact dependence of the shape 
on both strength and cutoff length is slightly dependent on back- 
reaction from modified expanding background, a detection of de- 
viation from the inverse-square law should be independent of the 
model for the cosmic expansion rate. 

The degeneracies discussed suggest caution when determining 
limits on deviation from the inverse- square law from galaxy clus- 
tering using sharp priors. At the same time, constraints on cosmo- 
logical parameters, such as Q.mo and /„, could be affected by mod- 



ified gravity physics. Full information from acoustic oscillations 
and evolution could disentangle some degeneracies. Furthermore, 
even if many imprints due to modified gravity are degenerate with 
other effects with respect to the shape, they are not with respect 
to normalization. While hot dark matter suppresses galaxy cluster- 
ing on small scales, enhanced gravity boosts perturbation growth 
on large scales. At the same time, in a cosmic expansion enhanced 
by enhanced gravity, perturbations are strongly suppressed and the 
standard growth rate is achieved only on large scales. 

We have discussed deviations from the inverse-square law in 
the standard framework of adiabatic initial conditions but it re- 
mains an open question whether a modified gravitation al forc e 
could make room for isocurvature initial conditions iPeeblesfcOoi) . 

Following a phenomenological approach, a Yukawa-like mod- 
ification to the potential can be introduced without being derived 
from some specific gravity theory. Although such an empirical 
analysis can be enough to investigate the effects of deviations from 
the inverse-square gravity law, some further theoretical considera- 
tions can be of interest. Variation of the gravitational constant can 
arise naturally in scalar-tensor theories. Whe reas a massive scalar 
field c an produce a Yukawa-like modification ( Z hvtnikov & Nesterl 
1 1994) . the existence of a massless scalar field coupled to the ten- 
sor field of Einstein gravity can implement a time varying G. 
The Jordan-Brans-Dicke theory is the original and simplest ex- 
tende d theory of gra vity which leads to time variations in G 
iBrans&Dicfall 19611) and has provided a suitable framework for 
inve stigation of the effects of a changing G on structure forma- 
tion <Liddle et alll998tlNagata et all200 2). The main effect on the 
matter power spectrum is connected to a shift to a higher redshift 
of the matter-radiation equality which on turn leads to a shift of the 
maximum of the matter power spectrum and to more small-scale 
power. 

A time-varying gravitational constant has an effect quite dis- 
tinct from that of a spatially varying gravitational constant, as con- 
sidered in this paper in the form of a Yukawa-like contribution to 
the potential. In fact, due to the lack of a characteristic scale dur- 
ing matter domination, a time-varying G will alter the growth rate 
but it will not change the shape of the spectrum, unlike a spatially- 
varying G, which naturally introduces a cutoff length scale. This 
main difference characterises potential parameter degeneracies. As 
we have seen, a Yukawa-like modification can act like hot dark mat- 
ter; on the other hand, the shift in the power spectrum due to a time- 
varying G has a potential degeneracy with increasing the number 
of massless species. Effects of either spatial or temporal variations 
in G could be better probed via their effects on the overall normal- 
ization of the power spectrum and its evolution, although biasing is 
a major source of uncertainty. 

Constraints on a Yukawa-like contribution to the potential 
have been found in the laboratory and analyses of planetary motion 
in the solar system can extend this to A ~ 1 AU ( Aldelberaer et al. 
2003). Such bounds on Aa can be very tight, but they do not ap- 
ply to scales beyond the solar system and give no information for 
A > 10 Mpc, i.e. the length scales we have been considering with 
regards to structure formation. On these scales, cosmology offers 
the only probe of modified gravity models. At present, we have seen 
that the near-degeneracy between modified gravity and changes in 
the cosmological parameters means that exotic gravitational mod- 
els remain hard to exclude. However, such degeneracies can be bro- 
ken with good enough data (e.g. Sawicki & Carroll 2005) and we 
can look forward to increasingly stringent tests of Einstein gravity 
on the largest scales. 
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